arXiv:1507.08937vl [cs.NE] 31Jul2015 


Efficient and robust calibration of the Heston option pricing 
model for American options using an improved Cuckoo 

Search Algorithm 

Stefan Haring Ronald Hochreiter 

Angnst 3, 2015 


Abstract 

In this paper an improved Cuckoo Search Algorithm is developed to allow for an efficient 
and robust calibration of the Heston option pricing model for American options. Calibration 
of stochastic volatility models like the Heston is significantly harder than classical option 
pricing models as more parameters have to be estimated. The difficult task of calibrating one 
of these models to American Put options data is the main objective of this paper. Numerical 
results are shown to substantiate the suitability of the chosen method to tackle this problem. 
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1 Introduction 

The classical textbook example for pricing options is the famous Black-Scholes model, see [1]. 
Since its creation in 1973, it has caused a dramatic increase in options trading because of its 
relatively simple usability. It provides closed-form solutions for European Call and Put options 
as well as American Call options. For American Put options, no closed-form solution exists 
because an optimal stopping time problem has to be solved. 

Over the years, criticism of the Black-Scholes model has arisen since it does not accurately 
capture the behaviour of option prices in the market. For example, the well-known volatility 
smile, which shows that deep in the money and out the money options have higher prices (and 
therefore higher implied volatility), is not taken into account in the Black-Scholes model since 
one of its main assumptions is constant volatility. 

To accurately reflect market behaviour, option pricing models in which volatility changes over 
time, so-called stochastic volatility models, have been created. In these models, the volatility 
is itself following a process that changes with time t. Some examples for these models are the 
Heston model [5], the CEV model [3] m and the Chen model [5]. 

Calibration of these stochastic volatility models is significantly harder than for the Black- 
Scholes model, since more parameters have to be estimated. The difficult task of calibrating one 
of these models to American Put options data is the main objective of this paper. We will focus 
on the Heston model as it is widely known and cited. In his seminal paper, Heston specified his 
stochastic volatility model for option prices as following: 

= ^dt + ^/V(t)dWi (1) 
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dV (t) = k(0 — V {t))dt + a^/V{t)dW 2 , (2) 

where Wi and W2 are Wiener processes that are correlated by dWi * dW2 = pdt. The 
parameters used for obtaining the variance process are the long run variance 6, the mean reversion 
rate k (the rate at which the volatility reverts to the long run variance) and the volatility of 
variance, a. 

Pricing options under the Heston model is not straight-forward when it comes to American 
option prices. The price of an American Call option is the same as for a European Call (because 
an early exercise is never optimal), however for American Put options no closed-form solution 
exists and therefore one has to use numerical approximations. There are many different methods 
in the literature for the approximations of these options such as the GARCH lattice approach 
as proposed by Ritchken and Trevor [6], tree-models presented by Beliaeva and Nawalkha [7] or 
Monte Carlo Simulation developed by Longstaff and Schwarz [8]. 

What is more important than choosing one of the many possible ways of how to price Amer¬ 
ican options under the Heston model is an efficient and robust way to calibrate the model to 
actual data. After all, using the model in practice for means of forecasting and pricing is the 
main objective. The goal of this paper is to fill the gap existing in research carried out on the 
Heston model. At the time of writing, no rigorously tested method of calibrating the model for 
American Call options has been published. 

The main reason is that while pricing options under this model is already a difficult task 
in itself, calibrating the model to a set of given option prices is even more challenging. In 
the standard model (keeping the interest rate hxed), there are five parameters that have to be 
calibrated: the aforementioned n, 9 and a as well as the starting value for the variance process 
Vq and the correlation between the stock price and the variance, p. 

What has been done, however, is calibrating the Heston model for European and Asian op¬ 
tions. Gilli and Schumann show the effectiveness in using Heuristics and Evolutionary Algorithms 
when calibrating the Heston Model for European options i, m- Collier gives an overview of 
algorithms that can be used for optimising a noisy objective function, in this case an Monte 
Carlo simulation, to calibrate the Heston model to Asian options data m- 

In this paper a naturally-inspired algorithm, namely Cuckoo Search, will be used to tackle 
the calibration problem. Naturally-inspired algorithms (also known as Evolutionary Algorithms) 
have been used on a lot of difficult optimization problems in Finance and have shown great 
promise (see for example na, m. m, usi)- See also relevant examples published in this 
journal including [16], [E], [E], [E], and [20] . 

The paper is structured as follows. Section gives a detailed overview over the problem and 
section will provide a brief introduction to the Cuckoo Search Algorithm. In section the 
chosen implementation for the algorithm will be explained while in Section numerical results 
will be presented. Section concludes this paper. 


2 Calibrating pricing models 

When calibrating a pricing model, the task is to find the parameters that minimize the differ¬ 
ence between the model price and the actual observed price. This is reflected in the following 
optimization problem: 


N 

min^(P™°‘^"' 

Z=1 
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(3) 
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Figure 1: Search space when varying volatility of volatility and rate of mean reversion 



volatility 


Figure 2: Search space when varying volatility and volatility of volatility 


where N is the total number of observed prices. The optimization problem can also be 
specified differently (mi), however the approach given in Equation (|^ is used as the objective 
function throughout this paper. 

The search space for the optimization problem is bounded continuous and not convex, which 
means standard optimization methods (e.g. gradient methods based on derivates) are not appli¬ 
cable. It is also characterized by a very high amount of local minima which complicates the task 
at hand even more. Figures andshow the results of the objective function given in Equation 
([^ when varying certain parameters. 

As can be seen from the figures, the optimization method must find the perfect balance 
between exploration (traverse the whole search space relatively quickly) and exploitation (search 
the near-neighbour regions of the best found solutions so far) to avoid getting trapped in the 
many local minima. 

In order for the algorithm to be a viable option, it has to fulfill the exploration-exploitation 
criteria as well as be numerically stable, meaning the error from the obtained solution and the 
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begin 

Objective function f(x), x = {xi ,..., xj)'^ 

Initialize population of n host nests Xi (i = 1, 2, ..., n) 
while stopping criterion is not met 

Get a cuckoo randomly by Levy flights 
Evaluate its fitness 
Choose a nest j among n randomly 
if (F, > F,) 

replace j by the new solution 

end 

A fraction (pa) of worse nests are abandoned and replaced by new ones 
Rank solutions and find the current best 

end while 

Postprocess results and visualization 

end 


Table 1: Cuckoo Search pseudo-code 


actual one has to be minimal. Another important criterion is computational efficiency. The 
algorithm has to converge to the solution as quickly and as stable as possible. 

The complexity of this task for the Heston pricing model for American options has so far 
prevented an effective and reliable method to arise. To this date, no algorithm has been proposed 
to solve this specific problem. In order to fulfill all of the above-mentioned criteria, a heuristic 
method called Cuckoo search will be used to calibrate the Heston option pricing model to a set 
of sample data. 

3 Cuckoo Search 

The Cuckoo search algorithm is a naturally-inspired optimization algorithm developed by Yang 
and Deb m- It is a so-called meta-heuristic algorithm meaning that it makes only few as¬ 
sumptions about the problem at hand and can be used effectively when there is only incomplete 
information and limited computational capacity given. This makes it an ideal candidate for the 
optimization problem outlined in section 

The heuristic used in this paper draws its inspiration from the behaviour of cuckoos that lay 
their eggs in nests of other birds that they come across during their flights. This so-called brood 
parasitism is translated into a very efficient optimization method. The cuckoo search algorithm 
follows three simple rules: 

1. Every cuckoo lays one egg in a random nest. 

2. The best eggs are carried over to the next generation. 

3. Cuckoo eggs are discovered by the host bird with a probability Pa G [0,1]. 

Here, the nests represent the amount of different solutions at each iteration of the optimization 
process. The eggs represent the solutions to the objective function and are obtained from the 
flight paths of the cuckoos. At each iteration there is a probability that a new solution gets 
discarded and replaced by a random result. The pseudo-code for the standard Cuckoo Search 
algorithm is given in Table 
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Figure 3: Levy flight in 2 dimensions with 1000 steps 


To model the flight paths that lead to new solutions, Levy flights are used (for an introduction 
to Levy flights see [22]) There is a plethora of research available on how the flight behaviour of 
insects, animals and even humans show the same characteristics as Levy Flights (see for example 

m. m)- 

A Levy flight is essentially a random walk where the step lengths are drawn from a Levy 
distribution, which is heavy-tailed. When the Levy flight is defined in a dimension greater than 
1, the steps are isotropic and random,meaning they are uniform in all orientations. 

The distinct feature of a Levy flight is made visible in Figure]^ Levy flights are characterized 
by large jumps followed by a series of small steps around the end of the jumps. This is exactly 
the behaviour needed for our optimization algorithm since it presents the perfect mix between 
exploration (the long jumps) and exploitation (the series of small steps). The long jumps also 
guarantee that the algorithm won’t get trapped in the many local minima present in the search 
space. 

A detailed explanation of the implementation of the Cuckoo search algorithm used in this 
paper is presented in the next section. 


4 Implementation 

The implementation of the Heston pricing model itself will not be discussed here since any 
method that fulfills the criteria of robustness and numerical stability as well as fast computation 
time is suitable. For this to achieve, the pricer should be implemented in C-|—h, Julia, Python 
or any similarly fast language. 

The Cuckoo Search algorithm used in this paper has been implemented as an R [25] package 
by the authors. Runtime tests using the package on standard optimization problems (Schwefel’s 
Function, Rastrigins Function, Easoms Function) were solved in sufficient time for the Finance 
application. 

In its essence, the implemented algorithm follows the steps laid out in the pseudo-code given 
in Table [^ It is modified in some parts to increase its effectiveness and convergence. 

The implemented method requires a small number of parameters because of its heuristic 
nature. The necessary input includes the number of nests (starting points and different solutions 
at every iteration), the number of iterations (or the error tolerance in the objective function). 
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Vvo 

0.2 

0.5 

0.3 

0.4 

a 

0.1 

0.1 

0.25 

0.25 

K 

3 

3 

2 

1 

0 

0.04 

0.25 

0.09 

0.16 

p 

-0.1 

-0.5 

-0.1 

-0.2 


Table 2: Parameter sets for the Heston pricing model 


upper and lower bounds for the parameters of the Heston pricing model as well as values Pmax, 
Pmin and Phase that are needed for the calculation of Pa ■ 

The part where the chosen implementation differs from the original algorithm is in the way Pa, 
the probability that the host bird discovers the cuckoos egg, is set. In the original paper, this is 
a fixed value. However Valian, Tavakoli, Mohanna and Haghi [55] show some modifications that 
improve accuracy and convergence of the cuckoo search algorithm. One of these modifications 
is that Pa is not fixed but a decreasing function between two set values, Pmax and Pmin- The 
function looks as follows: 


Pa{iter) = Pmax “ x ztcr, (4) 

where N is the total number of iterations and iter is the current iteration. This guarantees 
that the algorithm finds many new solutions in the beginning (the exploration is enhanced) and 
gradually moves on to hnd new solutions in the neighbourhood of points with an already good 
fitness (the exploitation is enhanced). Here we set the value for Pmax to 0.95 to have a lot of 
new solutions in the beginning of the optimization to further increase the exploration. The value 
for Pmin is Set to 0.05 so that at the end of the optimization the algorithm searches almost 
exclusively around the best found positions. 

The fitness function used to evaluate each solution in the search space calculates simply the 
sum of squared errors between the prices calculated by the Heston model and the observed (real) 
prices. It takes as inputs the parameters given by an egg in the Cuckoo Search algorithm and 
calculates the prices for all the asset value/strike price combinations. 

5 Numerical Results 

To test and calibrate the model, an artificial dataset has been used to validate the solutions. The 
parameters used in creating the dataset have been chosen to be similar to the dataset in Natural 
Computing Vol 4, Chapter 2 m- The model has been tested varying the spot price Sq as well 
as the strike price, K. In both occasions, the hxed price is set to 100 and the other price varies 
between 80 and 120 in steps of 2. 

The risk-free interest rate has been set to 5%, however it is very easy to adapt the model to 
allow for a varying interest rate as well. The maturities r for the options are q, | and 1 year. 
This results in a total amount of 84 prices that are used to calibrate the model. 

The parameter sets of the Heston pricing models are given in Table 

Lower and upper bounds for all 5 parameters based on logic and empirical findings have been 
introduced to restrict the search space only to plausible solutions. The boundaries are as shown 
in Table m 

For the correlation parameter, p, that captures the correlation between volatility and the 
stock price, the range has been set to be mostly negative according to several papers focused on 
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a 

n 

e 

p 

Lower Bound 

0.01 

0.05 

0.5 

0.01 

-0.8 

Upper Bound 

0.1 

0.3 

4 

0.1 

0.1 


Table 3: Lower and upper bounds for the parameters. 



iterations 


Figure 4: Convergence comparison for different nest sizes. 


this topic, e.g. Nandi m- 

For the choice of the number of nests used, multiple values have been tested. The algorithm 
has proved itself being very stable to changes in the number of nests, meaning values between 
20 and 50 have all shown fast convergence and high accuracy. 

From Figure]^ it is visible that for all choices of nests, convergence of the objective function 
is given and is relatively fast. This means there is no need to fit this parameter for the model 
and its selection can be made entirely to increase computational efficiency. The selection made 
for this paper is a nest count of 20. 

Furthermore tests on how many prices the algorithm needs as input for calibration have been 
conducted. Figure shows the convergence on a log-scale for inputs of 5,13 and 21 prices. It is 
clearly visible that the algorithm needs fewer iterations if the number of prices increases. There 
is, however, an increase in computation time for more prices. 

The graphs shown in Figure]^ show the convergence of the objective function as well as of all 
5 parameters calibrated by the Cuckoo Search algorithm. It can be seen, that there is a fast and 
stable convergence not only of the objective function but also of the parameters. The values of 
the objective function and the parameters after conducting the optimization process for each of 
the parameter sets are given in Table 


6 Conclusion 

In this paper an improved Cuckoo Search Algorithm has been developed to allow for an efficient 
and robust calibration of the Heston option pricing model for American options. Calibration of 
stochastic volatility models like the Heston is significantly harder than classical option pricing 
models as more parameters have to be estimated. The difficult task of calibrating one of these 
models to American Put options data has been solved using an implementation with the statisti- 
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iterations 


Figure 5: Convergence comparison for different number of prices 




Figure 6: Convergence of objective function and parameters. 


Objective function 


a 

K 

0 

P 

0.000001 

0.2 

0.1 

2.999 

0.04 

-0.09998 

- 

0.2 

0.1 

3 

0.04 

-0.1 

0.000001 

0.49999 

0.09999 

3.0012 

0.25 

-0.5 

- 

0.5 

0.1 

3 

0.025 

-0.5 

0.000001 

0.3 

0.250098 

2 

0.089978 

-0.1 

- 

0.3 

0.25 

2 

0.09 

-0.1 

0.000001 

0.39998 

0.25 

0.99875 

0.16 

-0.200125 

- 

0.4 

0.25 

1 

0.16 

-0.2 


Table 4: Objective function and parameter values after optimization. 












cal programming langauge R. The numerical results are very promising such that future research 
will focus on applying the algorithm to real-world Options data. 
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